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■ Abstract 
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, In this paper we try to estimate the lower critical dimension for replica 

symmetry breaking in spin glasses through the calculation of the additional 
free-energy required to create a domain wall between two different phases. 

ON . This mechanism alone would say that replica symmetry would be restored at 

■M ' the lower critical dimension D r = 2.5. 
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The calculation of free energy increase due to an interfaces is a well known 
method to obtain information about the lower critical dimension for spontaneous 
symmetry breaking. We perform the first analytic computation of this free energy 
increase in spin glasses and we use it to suggest the value of the lower critical 
dimension. 

Let us introduce the basic definitions. In the simplest case we can consider a 
system with two possible coexisting phases (A and B). We will study what happens 
in a finite system in dimensions D of size M d L with d — D — 1. Let us assume that 
in the d transverse direction we have periodic or free boundary conditions, while in 
the other direction (which we call t), we put the system in phase A at t — and 
in phase B at t = L. The free energy of the interface is the increase in free energy 
due to this choice of boundary conditions with respect to choosing the same phase 
at t = and t = L (for simplicity we suppose that the two choices A or B lead to 
the same value for the free energy). 

In many cases we have that the free energy increase SF(L, T) behaves for large 
M and L as: 

6F(L,M) = M d /L UJ (1) 

where u is independent from the dimension. There is then a critical dimension at 
which the free energy of the interface is finite: 

D c = u+1. (2) 

Heuristic arguments which sometimes can be made rigourous , tell us that at 
this dimension (the lowest critical dimension) the two phases mix in such a way that 
symmetry is restored. 

In most cases the value of u from mean field theory is the exact one and therefore 
we can calculate in this way the value of the lower critical dimension. The simplest 
examples are the ferromagnetic Ising model u = and the ferromagnetic Heisenberg 
model uj = 1. 

In this note we study the problem for spin glasses. For convenience we consider 
two replicas of the same system described by a total hamiltonian: 

H = H[a] + H[s] (3) 

where H is the Hamiltonian of a short range spin glass. 

In this case the order parameter is the local overlap = (<7jSj) between the two 
replicas. In mean field theory one finds that this overlap is constant in space, and 
all values in the interval [q m in, Qmax] are possible. A perturbative expansion around 
the solution with q constant in space can be done in sufficiently high dimensions, 
while infrared divergences appear in low dimension. 
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The aim of this note is to compute the free energy increase corresponding to 
imposing an expectation value of q equal to pi at t = and p 2 at t = L. We 
consider here only the case where both p's and non zero and are in the interval 

min j Qmax \ ■ 

We find that 

5F = V{\ Pl -p 2 \/Lfl\ (4) 

with V = M d L. As a consequence, the naive prediction of mean field theory for the 
lower critical dimension for spontaneous replica symmetry breaking is D c = 2.5. If a 
similar value would have been obtained for the interface free energy in zero magnetic 
field and with pi and p 2 of opposite sign, one would argue that the EA parameter 
should vanish at D = 2.5. 

We stress that these predictions are naive; corrections to the mean field theory 
are neglected. While in ferromagnets there is regime (low temperature) at which 
these corrections can be shown to be small, in the spin glass case the corrections to 
the mean field theory do not vanish even at zero temperature. We do not have here 
a regime in which we can show that the corrections to the free energy interface do 
not change qualitatively its L dependence. 

A simple testable prediction of our computation is that on a L D system with 
boundary condition q = qEA-, the expectation value of q in the centre of the box 
should go for large L as 

q (L) = q EA - const /LP ^ (5) 

Our paper is organized as follows. In section 1 we briefly outline the results of 
the mean field theory for short range spin glasses and we discuss its extension to 
the case of two constrained real replicas. In section 2 we show that equation (4) can 
be derived using rather general scaling arguments under the technical assumption 
that the overlap between the two constrained replicas varies linearly in space. A 
confirmation for that behaviour is found in section 3 where we find a variational 
approximation to the free energy increment SF. In the final section we draw some 
conclusions. 



1 The model 

The model we consider is the standard D-dimensional Edwards- Anderson spin-glass 
H on a square lattice in a finite volume V, which for simplicity will be taken as a 
box of size L » 1. This is defined by the Hamiltonian: 

H l S i] = ~ J ij S i S 3 ~hJ2$i (6) 

<ij> i 
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where with standard notations we have denoted by < i,j > the nearest neighbours 
on the lattice. The spins are Ising variables s» = ±1, and the couplings Jy are 
independent gaussian variables with zero average and variance Jfj = J 2 = 1. h is 
the magnetic field. 

In our discussion we will make extensive use of the results of the mean field 
theory (MFT) of the model near the critical temperature. Let us summarize here, 
without derivation, the main results. For a more complete exposition of the theory 



The study of the equilibrium properties of the model can be performed in the 
frame of replica method. The relevant order parameter is a space dependent over- 
lap matrix Qf 3 = {sfsf) which, analogously to the long range case, describes the 
statistics of overlaps between pure states. In MFT, where the system is treated in 
the saddle point approximation, one finds a (de Almeida-Thouless) line of second 
order phase transition to a glassy phase, which terminates for h = in T c — 1. 
In the vicinity of T c , the free energy as a functional of the order parameter admits 
a Landau expansion in which the original lattice is coarse grained, and one con- 
siders the order parameter averaged in small regions of the space V x centered in 
x, Qab{x) = j^-r J2 x &v x Qi b - The free energy functional in terms of Q a b(x) is then 
written as: 



where A is the Laplacian operator, r = T c — T = 1 — T, y = 2/3. The integration 
extends to the square box of size L, and 'Tr' denotes the trace in replica space. As 
usual, among all the quartic terms in Q which should be written in the expansion, 
we have only included the one responsible for the de Almeida-Thouless instability 
[0]. This leads to the phase transition into a replica symmetry breaking phase. 

The main advantage of this reduced model is that it allows for a complete analysis 
of the r.s.b.: the saddle point equations can be solved exactly above and below the 
transition temperature. In the low temperature phase, the solution to the saddle 
point equations is found in the framework of Parisi Ansatz ||, which in the present 
context consists in parametrizing the space dependent matrices Q a b(x) as a space 
dipendent functions q(x,u) with < u < 1. For free boundary conditions, the 
relevant saddle point is found to be constant in space, and the analysis become 
identical to that of the long range SK model |7[] . At the saddle point one finds: 



see e.g. @,§§§. 



2nF = Jd D x[±TrQ{x)AQ(x) + tTtQ 2 (x) + (l/3)TrQ 3 (x) 




(7) 
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with u = 3yq min , m = 3yq max . q min and q max are specified by the relations 

^yqlnin = h 

7~ qmaxiX ~2~qmax) (9) 

As in long range models, the appearence of r.s.b. imply the existence of many pure 
states with a non-trivial distribution of the mutual overlaps Q. The theory predicts 
that these overlaps are constant in space. 

Let us now enter into the core of our discussion. We want to introduce boundary 
conditions in the model to force spatial dishomogeneity of the order parameter. In 
analogy with what is done in ordered systems we would like to put the system in 
two different equilibrium states at the two boundaries along a given direction. We 
observe that this can not be done by imposing boundary values to the function 
q(x,u). Any variation with respect to the form would take us outside of the 
saddle point, where the free energy functional has no physical meaning. We shall 
follow instead a procedure introduced in a previous paper |T0| (referred as I in the 
following), to study long range spin glasses off-equilibrium. 

Consider two identical (same Jy) replicas of the system, which undergo different 
thermal histories (for equal temperature and magnetic field). We constrain the 
overlaps between these two real replicas (RR in the following) on the two boundaries 
along one direction and leave free boundary conditions on all others. In this case 
the Saddle Point overlap will be constant in all but the privileged direction, and 
we will have to solve an effective one dimensional problem. Denoting d\ and 82 the 
boundaries on which we impose the non trivial constraint, we can write the partition 
function for the doubled system as 

Z c = £ exp {-/3H[s] - /3H[<t]} II TtTT £ " ^ 

■ nMipr ( 10 ) 

The object of interest will be the free energy difference SF between this situation 
and the unconstrained case where the delta functions are not present in the partition 
sum. We will get an estimate of SF in mean field theory. This will enable us to 
estimate the probability of fluctuations of the overlap of amplitude \p\ — P2I over 
a scale L through Prob(\pi —p2\,L) ~ exp(— f35F). If this probability remains 
finite for large L, this kind of fluctuations destroy replica symmetry breaking. In 
high enough dimension we will find a free energy difference divergent with L. The 
critical dimension D c will be then identified as the dimension at which this free 
energy difference become finite. 
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In the replica treatement of the problem, one has to replicate both the s and the 
a spins [|K], |TT| . Thus three space dependent n x n matrices will appear: QH(x) 
describing the overlaps among s spins, QH(x) describing the overlaps among a 
spins and Qali x ) describing the overlaps between s and a spins. For symmetry 
reasons Q 12 b [x) = Q b l(x); the diagonal elements Q aa are as usual taken to be zero 
by convention. We will choose in the following saddle points for which Q l J b {x) = 
Qab(x) = Qab(x) and Qab(x) = Qll(x) = P a b(x) . The constraint introduced in the 
partition function reflects itself in the order parameter through the fixing of the 
values of the diagonal elements of Q^l on the border: 



All other elements of the matrices are to be determined from the variational equa- 
tions for the free energy. In dealing with the matrices Q r a l r,s = 1, 2 a,b = 1, n 
it is useful to introduce new indices a = (r, a), (3 = (s,b) and a In x In matrix 
Qa/3 = Qat that contains all the three matrices. In term of this new matrix, the free 
energy functional is formally identical to that for a single real replica in term of the 
usual Q ab . The difference lies in the variational procedure, where one has to keep 
into account the constraint (|TTp . 

Near the critical temperature, the free energy admits a Landau like expansion 
(see 0),with the single replica matrix Q substituted by Q. The saddle points equa- 
tions in terms of the matrices Q and P are: 



As usual, to solve these equations we need an ansatz that will eventually allow us 
to do the analytical continuation to n — > 0. As in I we assume both matrices Q 
and P to be hierarchical matrices. Namely, we parametrize Q a b{x) by a function of 
u G [0,1] q(x,u) and the matrix P a b(x) by a diagonal element P aa (x) = p(x) and 
a function p(x,u). Physically p(x) represents the space dependent overlap between 
two RR constrained on the boundary. A tentative discussion of the physical meaning 
of the functions q(x, u) and p(x, u) can be found in I. 

Clearly, for the chosen boundary conditions, the various parameters will be con- 
stant in all but one direction. Labelling the coordinate along this direction by t we 
find that q p and p will depend on space only through t. In this way, denoting the 
integrals of the kind Jq 1 du g(t, u) as (g), the saddle point equations become: 



(11) 



AQ ab 
AP ab 



2rQ ab + (Q 2 ) ab + (P 2 ) ab + yQl b + h 2 
2rP ab + 2(QP) ab + yP^ b + h 2 



(12) 



dt 2 



q{t,u) 



2(r-(q))q(t,u) + 2(p-(p))p(t,u) 
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at 



2 p(t, u) 



d 2 

Of 



+ / dv [q(t,u) -q{t,v)} 2 + / dv \p(t, u) - pit, v)} 2 + yq\t, u) + h 2 
Jo Jo 

= 2{r-{q))p{t,u)+2{p-{p))q{t,u) (13) 
dv [q(t, u) - q(t, v)] \p(t, u) - p(t, v)] + yp 3 (t, u) + h 2 

2rp{t) -2(qp) + yp{tf + h 2 
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If pi = p 2 and q min <P\< q ma x one can find a solution without spatial dependence 
at all, with p = p\. Then the problem reduces to the one discussed in I, There 
we showed, on very general grounds, that for any p\ in this interval, there exist a 
solution to the saddle point equation which has the same free-energy density of the 
unconstrained system (5F = 0). The results of I in the present context are: 



q{x,u) 



q(u) 



p(x, u) = p(u) 



Qmin 


U < Mo/2 


2x/(3y) 


Mo/2 < u < Up/2 


< p 


u p /2 < u < Up 


x/3y 


Up < U < U\ 


qmax 


U\ < u < 1. 


qmin 


u < uq/2 


2x/(3y) 


uq/2 < u < Up/2 


[ V 


Up/ 2 < u < 1 



(14) 



p(x) = p = p\ 



(15) 
(16) 



The parameters q m i n , q ma x, u and u\ are those of the unconstrained solution (|3|) and 
Up = 3ypi is the point where the function fl8|) is equal to p\. This solution reflects 
the fact that the multiplicity of states does not give an extensive contribution to 
the free energy. The space of equilibrium states of the two copies constrained at an 
overlap p\ is simply a restriction of the cartesian product of the equilibrium states 
of two free system. The freedom in the choice of p±, which is a zero mode of the free 
energy, is the spin glass analog to the Goldstone zero mode found in ordered models 
with a spontaneously broken continuous symmetry. It is now clear that if we impose 
Pi 7^ P2, but both in the interval [q m i n , q m ax}, we are choosing at the boundaries two 
of the possible overlaps admitted by the free problem. The additional free energy 
will have to go to zero in the limit when the boundaries become very far from each 
other. Our aim is to know with what power it goes to zero. 
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2 A dimensional argument. 



To study the situation pi ^ p 2 we could perturb around the solution with p\ = p 2 . 
Above the critical dimension the laplacian term in (|7|) can be treated as a small 
perturbation. The relevant parameter of the expansion turns out to be \p 2 —pi\/L 
which is always arbitrary small. 

It is reasonable to think that the solution to the saddle point equations would 
be in this case similar in form to eq. (|T5|) , but with p = p(t) a function interpolating 
between pi at t — and p 2 at t = L. On physical grounds we expect that after av- 
eraging over the quenched disorder, p(t) interpolates linearly between the boundary 
values, namely 

p{t)=p 1 (l-t/L)+p 2 t/L. (17) 

This assumption will enable us to determine the lower critical dimension by mere di- 
mensional analysis. Our hypothesis will be validated a posteriori in the next section, 
where we will derive the linearity in the context of an approximate maximization of 
the free energy. 

To first order in the perturbative expansion in pi — p 2 , the variation of the free- 
energy is n5F = J dt TrQAQ computed at the unperturbed saddle point. It is easy 
to see that this variation is zero for the saddle point (|T5*|): 

I /dt ^Q q/3 AQ q/3 (18) 

= J dt (^^j [1 - / du9[u - u /2)6(u - u) - j du6(u - u Q /2)\ = 0. 

A non zero 5F at this level would have implied, by simple dimensional analysis, 
5F ~ L D ~ 2 , i.e. the same result found for ordered systems with a continuous 
symmetry. The vanishing of this term means that D c is higher than 2. One could 
say that fluctuations due to the zero mode in spin glasses cost less and become 
therefore important earlier than the usual Goldstone modes in ordered systems. 

Substituting the expresion (pl|) into the saddle points equations ([H|) and denot- 
ing Up = 3yp, x — %(dp/dt) 2 we obtain 

d 2 g{t,u) 1 

— — — = x[o{u-u p ) - -5(u-u p /2)\ 

^('•") = - xS (u-« r /2). (19) 



dt 2 



We find in this way that the expresion ([15]) satisfies the saddle point conditions for 
all u except u p /2 and u p . It is reasonable to suppose that the effect of the Laplacian 
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term in the free-energy will result in a (small) smoothing of the functions q and p 
around u p /2 and u p . Thus we suppose that the functions will have variations of a 
given order e in regions of order e' around u p /2 and u p . e and e', as well as u p will be 
in general function of t. Using the monotonicity of the functional relation between 
p and t, we will consider all parameters as functions of p. 

We now show that under this hypothesis the following remarkable facts happen: 

• 8F, defined as F(j>2,Pi)—F(pi,pi) does not depend neither on the temperature 
r nor on the magnetic field h. 

• 5F must behave as L D x (5/4) , in order that the linear interpolation ([H]) maxi- 
mizes the free energy. 

Let us write the variational matrix Q as Q = + <5Q where solves the 
variational problem in the absence of the perturbation. 5F can be written as 

- 2n5F = [dx\~ Tr(25QAQ (0) + 5QA5Q) (20) 



2 

+ tTt5Q 2 + TrQ (0) <5Q 2 + ^TMQ 3 (21) 



a/3 



(22) 



Because of the unperturbed SPE the only terms linear in <5Q that can appear arise 
from the Laplacian term. For a given t (or equivalently p), one can evaluate the 
variation of the free energy density 5f just dropping the integration over x in (|22|). 

In this free-energy we want to study the dependence on r and h before optimizing 
with respect to <5Q. As one can see from (|T5|), the function p does not depend on 
r and the function q depends on r only in the region u > u max where q{u) = q max - 
In the same way one observes that both functions depend on h only in the region 
u < u min where they take the value q m in- To study the dependence of 5F on 
the temperature and the magnetic field we can use q max and q min as independent 
variables instead of using r and h which give rise to simpler algebra. The structure 
of 5Q we choose implies that the only terms which can depend on q ma x are: 

• rTr5Q 2 , trough its r dependence (r = g ma!r (l - 3yq max /2)), and 

• TrQ(°)<5Q 2 

The other terms do not depend on q max because for the a and (3 such that Q^°j = q max 
one has <5Q a/ 3 = 0. So we find that 

- 2n^- = (1 - 3yq max )Ti5Q 2 + £ 8Q 2 af) (23) 

a P\H)yk=<lmax 
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It is easy to show, using the algebra of ultrametric matrices, that the second adden- 
dum on the r.h.s. of fl23| ) is exactly equal to —(1 — 3yq max )TrSQ, 2 and consequently 
Sf is independent of r. Along the same lines one shows that the only possible de- 
pendence on q min is in the term TrQ(°)<5Q 2 , but its derivative w.r.t. q min is equal 
to zero. Therefore the free-energy density variation can only depend on p. Let us 
analize the dimensions of each term in Sf in terms of p, e, e' and x- We remind that 
while p and \ are fixed parameter in the problem e and e' are to be determined by 
sadlle point equations. In the considerations which follow we can safely assume, that 
dimensionally e' ~ e. Keeping the two quantities different would only complicate 
the formulae, but not the scaling of the free energy. 

We find for the dimensions of the different terms in fl2"2"|): 

Tr25QAQ (0) 

rTr<5Q 2 + TrQ (0) 5Q 2 

Tr5Q 3 ~ e 5 (24) 

a/3 
a/3 

a/3 

If we now rescale: 

e — > pr] (25) 
X ^ (26) 

Sf becomes an homogeneous function of order 5 in p. If the linear form p(t) has to 
be a maximum of the free-energy, the free energy density must be independent of 
p. This can be seen from the fact if one optimizes for fixed p(t) one finds in general 
Sf = Sf(p, (dp/dt)). The further minimization with respect to p(t) leads to Euler- 
Lagrange equations that give zero 'acceleration' d 2 p/dt 2 only if Sf is independent of 
p. So the p 5 dependence must be compensated at the saddle point by the dependence 
on cf) = x/P A - m this way, writing in all generality for the saddle point Sf = p 5 g(<f)) 
one must have g((f>) ~ c// 5 / 4 ) = X /P 5 - Using x = 3y(P2 —Pi)/L for the form (|17D , 
we find that the total free-energy variation scales as: 

SF ~ L D X m) ~ \p 2 -pi\ {5/2) L D - 5 / 2 . (27) 
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Equation ( P7| ) is the main result of this paper. It tells that in the context of mean 
field theory, the lower critical dimension at which SF become finite for finite \p2 —pi\ 
is D c = 2.5. Moreover for D > 2.5 one can expect fluctuations in space of the order 
parameter to scale with the distance as \q(x) — q(y)\ ~ \x — y\ l ~ 2D ^. It is worth 
noticing that the value for the critical dimension we get is fully compatible with the 
one found by Bray and Moore for the glassy transition in absence of field [|T2J with 
totally different methods. The same and other authors [[13], Q claimed that in any 
finite dimension the spin glass transition is destroyed by the presence of a magnetic 
field. Our findings disagree with this claim, as we find Sf to be independent on the 
magnetic field. 



3 A variational approximation. 



Let us now turn to an explicit computation of the free-energy density increment 
through a variational approach. Instead of solving the full SPE ([13]) we will here 
propose an explicit parametrization of the small variation to the form (15") in the 
neighbourhood of the points u p /2 and u p , and we will maximize the free energy with 
respect to the parameters of this variation. We expect that the numerical value so 
obtained for 5F is a lower bound to the real value. Furthermore we will show that it 
scales as discussed in the previous section. Thus this section proves that the leading 
behaviour previously obtained does not accidentally cancel. Moreover here we will 
not need to assume the linear form (|17D for p(t): this will be found as the optimum 
of the free energy. 

We choose to smooth the singularities around u p /2 and u v (cf. (|l^)) by interpo- 
lating with an arc of parabola the step-wise linear behaviour of the functions q and 
p in the surroundings of these points. Our variational functions will then be: 



q(u) 



p(u) 



qmin 




U < Uq 


2«/(3y) 




Uq < U < U\ 


p — a(u 


-u 2 f 


U\ < U < U 2 


p 




U,2 < U < M3 


p + a'(u 




U3 < U < W4 


u/(3y) 




U4 < U < U 5 


qmax 




U > U 5 


qmin 




U < Uq 


2u/(3y) 




Uq < U < U\ 


p — a(u 


-u 2 f 


Ui < U < U2 


p 




U > U 2 



[28) 



(29) 
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with 

u = 3yq min /2 u 1 = 3yp/2- S/2 u 2 = 3yp/2 + S/2 

u 3 = 3yp - 5'/ 2 u A = 3yp + S'/2 u 5 = 3yq max (30) 

a=l/(3y5) a' = 1/(6^'). 

The reader should not be confused by the notation at this point, although q min 
and q max are the same as in the presious sections and of (0), we changed here the 
notation for the points {% = 0, 5). The only variational parameters that appear 
in the the free energy functional are S and 8' . Denoting as in the previous section 

X = 3y f ^ , the free energy density increment as a function of S and S' takes the 
form: 



315 5 S A p xl $ 3 5' 2 S 2 p8' 2 pS' 4 118' 5 

102060?/ 3 324y 2 + 9y + lSy 9720y 3 + 6A8y 2 518V 3265920?/ 3 



which has to be maximized with respect to 5 and 5'. Equation (pTf) is consistent 
with the scaling established in the previous section. The change of variables: 

„a + b 

o = yp^— 

5' = yp b -^ (32) 

x = y 2 P 4 4> 



gives us 



Sf = y 2 p 5 (29 a 4 - 10080 a 2 b - 252 a 3 b 

J y F 52254720 V 



+ 142 a 2 6 2 - 84 a 6 3 - life 4 + 2903040 0). (33) 

It is apparent form the variational equations for a and b that a should be of or- 
der while b ~ 1 / 4 . To lowest order in (ft the solution is a = — 6 2 /240, 
b = 12v / 2(7/ll) < - 1//4 ' ) 1//4 which gives for the free-energy density 

Sf = 0.673659 y\ 5/4 . (34) 

This result confirms in a specific example the behaviour in x 5 / 4 of Sf which is the 
only one compatible with p(t) linear in t. The total free-energy of the interface, 
SF = SfL D is proportional to L D ~ 5 / 2 \pi — P2\ 5 ^ 2 , confirming the analysis of the 
previous section. Let us observe here the proportionality of SF to y 2 , the coefficient 
of the quatic term in (^). 
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4 Conclusions 



In this paper we have shown that it is possible to estimate the cost in free energy 
of a domain wall between two different phases in spin glasses. This information was 
used to indicate a possible value for the lower critical dimension D c = 2.5 for replica 
symmetry breaking. 

We are aware of the several criticisms that could be raised against this indication. 
In ordered system MFT gives a reliable estimate of D c because the fluctuation of the 
order parameter are negligible at very low temperature. In the spin glass case there is 
not such evidence. Moreover, in ordinary systems, this kind of analysis is confirmed 
by the behaviour of the perturbation expansion. In O(N) models, for example, the 
free propagator G(k) ~ k~ 2 in the ordered phase and therefore the fluctuations of 
the order parameter diverge in D = 2. In the spin glass case, there is a whole family 
of propagators, the most divergent of which, in presence of magnetic field, has the 
behaviour G{k) ~ k~ 3 [15]. If this behaviour is not modified by renormalization it 
would imply a lower critical dimension D c = 3. 

At present we do not know if any of this criticisms will be substantiated and the 
actual critical dimension in spin glasses is larger then 2.5. But if this happens, the 
replica symmetry restoration mechanism must be different from the simple 'instan- 
tonic' one that we have proposed in this paper. 

The present state of the art in numerical simulations of spin glasses indicates 
that in dimension two there is no transition |16| while in dimension four there is 
full replica symmetry breaking ||17|| . Unfortunatly the situation is far from clear in 
dimension three. Recent simulations by one of us |18| obtained from large lattices 
were compatible with a finite temperature phase transition but also with an essential 
singularities at T = 0, which would indicate D c = 3. 

All this calls for further research. One can expect that simulations on larger 
lattices will eventually resolve the problem of the existence of the transition in 
three dimensions. On the more analytical side a crucial problem to attack is the 
renormalization of the k~ 3 behaviour of the most singular propagator in perturbation 
theory. 
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